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Abstract. Here I summarize constraints on the nature of the dark matter halos of field galaxies that have been obtained from 
the most recent investigations of (i) weak galaxy-galaxy lensing and (ii) the dynamics of satellite galaxies in orbit about large 
host galaxies. Both of these techniques are statistical in their nature (i.e., large samples of galaxies are required to obtain 
a "signal"), but since they have inherently different selection biases and systematic errors, they are quite complementary 
to each other. Results of work over the last several years on weak lensing and satellite dynamics is revealing a remarkably 
consistent picture regarding the dark matter halos of bright field galaxies (L > V). The halos extend to large physical radii 
(> 150 h~ x kpc) and are flattened in projection on the sky, there is a marked difference in the depths of the potential wells of 
early-type galaxies and late-type galaxies, and the velocity dispersion profiles of the halos, <j v (r p ), decrease at large projected 
radii. All of these are expected to hold true in a cold dark matter universe and, while neither technique can address the possible 
small-scale 5 h^ 1 kpc) conflicts between cold dark matter and observed galaxies, on scales £ 50 h kpc both techniques 
yield results that are consistent with each other, and with the predictions of cold dark matter. 



INTRODUCTION 

The existence of dark matter halos surrounding large, bright galaxies is well established (e.g., (J], 01, @], 0] and 
references therein), and in the standard cold dark matter (CDM) paradigm, the halos of large field galaxies are expected 
to extend to virial radii of ~ 100/i~' kpc to ~ 200/1"' kpc and have masses of ~ \Q n hT l M & (e.g., 0|, @], 0]). Until 
very recently, however, direct observational constraints on the nature of the dark matter halos of field galaxies have 
not been especially strong. In particular, it has been challenging to address the question as to whether the halos of 
observed galaxies are consistent with the halos that one would expect in a CDM universe. 

The lack of a Keplerian fall-off in the rotation curves of the disks of most spiral galaxies (e.g., |2[) indicates 
that the dark matter halos extend far beyond the visible radii of the galaxies. Therefore, in order to place constraints 
on the total mass distribution, it is necessary to use tracers of the halo potential that exist at large projected radii 
(> 100/i~' kpc). Two such tracers of the large-scale potential are satellite galaxies that are in orbit about isolated 
host galaxies, and photons emitted by distant galaxies that, on their way to the observer, happen to pass through the 
potential wells of more nearby galaxies at small impact parameters (i.e., gravitational lensing). "Strong" gravitational 
lensing, in which multiple and highly-distorted images of a source occur, is a rare phenomenon because it requires 
nearly perfect alignment of the lens and source galaxy (e.g., |8|]). "Weak" gravitational lensing, in which multiple 
images and significant image distortion do not occur, is, however, commonplace in the universe (e.g., |9], 1 10], 1 11], 
fl2Tn . and it is on this extremely mild regime of gravitational lensing that I will focus for this discussion. 

Weak lensing of background galaxies by foreground galaxies ("galaxy-galaxy" lensing) and the motions of satellite 
galaxies about host galaxies are phenomena that can only lead to constraints on halo potentials through ensemble 
averages over statistically large samples. That is, for any given foreground galaxy, the distortion that it induces in 
the images of background galaxies due to weak lensing is so small that the signal cannot be detected convincingly 
for any one foreground lens galaxy. Similarly, isolated host galaxies are typically found to have 1 to 2 satellite 
galaxies on average and, so, the potential of any one host galaxy cannot be constrained at all well by the motions 
of its own satellites. Both galaxy-galaxy lensing and satellite dynamics, therefore, lead to statistical constraints on 
the halo population as a whole and by their very nature they require large samples of galaxies in order to obtain such 
constraints. Until several years ago, galaxy-galaxy lensing and satellite dynamics were both tantalizingly close (or 
frustratingly close, depending on one's point of view) to being able to fulfill their theoretical promise to map out the 
gravitational potentials of the halos of field galaxies. With the advent of routine availability of wide-field imagers 
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FIGURE 1. Schematic representation of multiple lenses along the line of sight to a given source galaxy. 



and the completion (or near completion) of large redshift surveys, however, both galaxy-galaxy lensing and satellite 
dynamics are now yielding sufficiently strong constraints on the dark matter halos of galaxies that the observations 
can be used to test the theoretical predictions (i.e., CDM) at a substantive level. 

There are distinct advantages and disadvantages of galaxy-galaxy lensing versus satellite dynamics when it comes 
to constraining halo potentials. A clear advantage of galaxy-galaxy lensing is that it can be applied to all foreground 
galaxies and, since gravitational lensing is affected only by the total mass along the line of sight and not its dynamical 
state, the halos of the foreground lens galaxies need not be virialized. A complicating factor in galaxy-galaxy lensing 
is that it is not correct to assume that each background galaxy has been lensed solely by one foreground galaxy 
(e.g., fPjll . fl4ll . (3)- Instead, photons emitted by the distant galaxies are deflected by all mass along the line of 
sight, including individual galaxies, groups, and clusters (e.g., Figure^. That is, galaxy-galaxy lensing is inherently a 
multiple-deflection problem and care must be taken when using observations of galaxy-galaxy lensing to constrain the 
halos of a given subset of lens galaxies (i.e., the halos of early-type galaxies versus late-type galaxies, or the halos of 
high-luminosity galaxies versus low-luminosity galaxies). Therefore, a computation of the weak lensing signal about 
the white lenses in Figure ^ above is not identical to a measurement of weak lensing signal produced by the white 
lenses since the black lenses also contribute to the net shape of the final image. That is not to say that galaxy-galaxy 
lensing cannot be used to probe the potentials of halos surrounding lenses of differing types; it most certainly can, but 
the presence of multiple deflections in the data must be taken into account when modeling the observed signal. In the 
case of relatively shallow data (zi ens ~ 0.15), most sources will have been lensed by only one foreground galaxy (e.g., 
1 16]), but in deep data sets (zi en s ~ 0.5) most source gala xies will have been lensed at a significant and comparable 
level by two or more foreground galaxies (e.g., 1 13], 1 15]). A further disadvantage of galaxy-galaxy lensing is that 
the signal is very small (systematic image distortions of < 1% or in the image ellipticities), so the images of millions 
of background galaxies must be obtained and, in general, be meticulously corrected for the presence of anisotropic, 
spatially-varying point spread functions. Finally, it is possible that Newtonian tidal distortions of genuine satellites 
of the lens galaxies could masquerade as a weak lensing signal. Happily, such distortions appear to be at most a very 
small contributor to the observed weak lensing signal (e.g., 1 13], 11711 . 11811 . I19I1 ). 

An advantage to using dynamics of satellite galaxies to probe the potentials of the halos of isolated host galaxies 
is that, unlike deep weak lensing data, the only important potential well in the problem is that of the host galaxy. 
In principle, this is a "cleaner", more straightforward probe of the halo potential which is intentionally restricted to 
the physical scales that one would expect to characterize the halos of individual large galaxies. However, there are 
a number of arguments against the use of satellite dynamics to probe the mass distributions of host galaxies: (1) 
satellites must be found at large projected radii in order to probe the halo potential on the very largest scales, (2) 
noise is introduced by the presence of "interlopers" (i.e., galaxies that are selected as satellites but which are, in fact, 
not associated dynamically with the host galaxy), and (3) the relaxation times of these systems are large compared 
to the age of the universe. The first argument is much less compelling now than it was in the past simply because of 
the availability of large redshift surveys (i.e., the data bases are now sufficiently big that although it is rare to find 
satellites at, say, a projected radius of 500 h kpc, the large number of redshifts that are now available makes it 
possible to compile statistically significant samples). The second argument has also become much less compelling 
with the realization that it is straightforward to account for the effects of interloper galaxies on the determination of 
the velocity dispersion (see below). The third argument can still be compelling, since it makes little sense to apply a 
virial-type mass estimator to systems which are not relaxed. However, an assumption of virialization is not necessary 
a priori, and the use of secondary infall models can be used to bypass this assumption (e.g., 1 20]). 

For the sake of a certain amount of brevity and, at the very least, an attempt at providing some level of coherent 
argument, I will focus here on only the most recent results that are directly relevant to the halos of field galaxies and 
which have been obtained from four large surveys: the COMBO-17 Survey, the Red-Sequence Cluster Survey (RCS), 



the Sloan Digital Sky Survey (SDSS), and the Two Degree Field Galaxy Redshift Survey (2dFGRS). Even with this 
restriction, it is simply not possible to discuss all of the most recent results from these studies in great detail, and the 
reader should consult the source literature for further information. Finally, I should hasten to add that all errors or 
omissions in this article in regards to my colleagues' work are entirely unintentional and entirely my own fault. I can 
only hope that my colleagues will be kind enough to forgive me. 



THE SURVEYS 
COMBO-17: Galaxy-Galaxy Lensing 

The acronym COMBO-17 stands for "Classifying Objects by Medium-Band Observations in 17 filters" l2lll . l22ll . 
ll23ll . The COMBO-17 survey consists of high-quality imaging data with the ability to obtain both rest frame colors 
and accurate photometric redshifts (5z/(l + z) < 0.01 for R < 21, 8z/(l +z) ~ 0.02 for R — 22, and <5z p hot < 0.1 
for R < 24). The survey consists of 5 fields, including an extended region in the location of the Chandra Deep 
Field South (CDFS). The observations were carried out using the Wide Field Imager at the 2.2-m MPG/ESO 
telescope. The field of view of the camera is 34' x 33' and the 17-band filter set covers a wavelength range of 
350 nm < A b s 5! 930 nm. The latter allows for a rough determination of the spectral energy distributions of the 
objects, which in turn leads to both reliable classification of the objects into galaxies, quasars, and stars, as well as 
the ability to determine accurate photometric redshifts. A catalog containing astrometry, photometry in all 17 bands, 
object classification, and photometric redshifts for the 63,501 objects in the extended CDFS is publicly-available 
lEl (see |http: //cdsweb . u-strasbq . f r/cqi-bin/qcat? J/A+A/421/913> . The COMBO-17 results 
that will be discussed here consist of efforts to use galaxy-galaxy lensing to study dark matter halos. The data set is 
particularly well-suited to this task because of the reliability with which background galaxies (i.e., lensed sources) can 
be separated from foreground galaxies (i.e., the lenses). Note, too, that although the full COMBO-17 survey covers 5 
fields, the results shown here come from only 3 of the fields (a field centered on the cluster A901, the CDFS field, and 
a random field 12410 . 



RCS: Galaxy-Galaxy Lensing 

The RCS (|25], ll26l0 is a somewhat shallow (5(J point source detection limits of Rc ~ 24.8 and z' ~ 23.6), wide 
field (~ 90 sq. deg.) imaging survey that was designed primarily to search for galaxy clusters out to redshifts of z ~ 1 .4. 
The images for the complete survey were obtained with the CFHT and CTIO 4— m telescopes using mosaic cameras, 
and consist of 22 widely-separated patches of ~ 2.1° x 2.3°. The RCS results that will be discussed here consist of 
galaxy-galaxy lensing studies and were obtained from ~ 42 sq. deg. of northern RCS data. Without spectroscopic 
or photometric redshift information, the RCS galaxy-galaxy results had to be obtained from a rough separation of 
lenses and sources that was based upon apparent magnitude cuts (i.e., galaxies with "faint" apparent magnitudes are 
on average background objects while galaxies with "bright" apparent magnitudes are on average foreground objects). 
Although the foreground-background distinction between a given pair of galaxies in the RCS data is by no means 
as secure as in the COMBO-17 data, the RCS is nevertheless a superb data set for galaxy-galaxy studies simply 
because of the area covered (~ 45 times larger than COMBO-17 for the weak lensing work). Given that weak lensing 
is primarily a statistical game, this is a good example of how well the galaxy-galaxy lensing signal can be detected 
and also used to constrain the nature of dark matter halos given only minimal distance information and a tremendous 
number of candidate lenses and sources. 



SDSS: Galaxy-Galaxy Lensing & Satellite Dynamics 

The SDSS is a combined photometric and spectroscopic survey that will ultimately map roughly one quarter of the 
sky above I ~ 30° and provide redshifts of ~ 10 6 galaxies and ~ 10 5 quasars with r' < 17.8. The SDSS is a fully- 
digital survey and makes use of 5 broad optical bands (u', g', r 1 , i', z') for photometry. The data for the SDSS are being 
acquired at the Apache Point Observatory in Sunspot, New Mexico using a 2.5-m telescope, as well as three, smaller 
subsidiary telescopes for the purposes of photometric calibration, monitoring of the seeing, and scanning for clouds. 



The rms galaxy redshift errors are ~ 20 km sec~'to ~ 30 km sec~'(e.g., ¥Fh . A technical summary of the 

SDSS can be found in York et al. l29Tl . information about the main galaxy sample is given by Strauss et al. 1 30], and 
information about the photometric system and photometric calibration is given by Fukugita et al. l3lll . Hogg et al. l32ll . 
and Smith et al. [33]. All of the SDSS data, including astrometry, photometry, redshifts, and spectra, are available via 
the SDSS website (http : / / www . sdss . orgi using structured queries that can search and combine the individual 
data bases. The third SDSS data release occurred on September 27, 2004 and includes spectra of 374,767 galaxies, 
spectra of 51,027 quasars, and photometry of 141 million unique objects. The SDSS results that will be discussed here 
consist of both galaxy-galaxy lensing studies and studies of the satellites of large, isolated galaxies. 



2dFGRS: Satellite Dynamics 

The 2dFGRS is a spectroscopic survey in which the target objects were selected in the bj band from the Automated 
Plate Measuring (APM) galaxy survey ( ll34Tl . O) and extensions to the original survey. A detailed discussion of 
the survey and the data base is given by Colless et al. |36]. The observations, which are now complete, were carried 
out at the Anglo-Australian Telescope using the Two Degree Field (2dF) multifiber spectrograph. The final data 
release occurred on June 30, 2003 |37] and includes reliable redshifts of 221,414 galaxies with extinction corrected 
magnitudes of bj > 19.45, covering an area over ~ 1500 square degrees. Galaxies with reliable redshifts have 
an rms uncertainty of 85 km sec -1 !^ . All data, including spectroscopic catalogs (245,591 objects), photometric 
catalogs (382,323 objects), and FITS files containing the spectra, are publicly-available from the 2dFGRS website 
Qhttp : //msowww . anu . edu/au/2dFGRS| l. The 2dFGRS data base is fully-searchable via structured queries, 
and on-line documentation is available on the 2dFGRS website. The photometric transformation from the SDSS band 
passes to bj is 

bj = g l + Q.\55 + Q.\52{g l -r l ) (1) 

(e.g., 1 38]). The 2dFGRS results that will be discussed here consist of investigations into the nature of the dark matter 
halos of large, isolated galaxies that are orbited by one or more satellite galaxies. 



PROBING HALO POTENTIALS WITH WEAK LENSING 

General Relativity tells us that any mass will cause a curvature of spacetime in its vicinity. Therefore, any mass 
located along the line of sight to a distant luminous object will act as a gravitational lens by deflecting light rays 
emanating from the object as they propagate through the universe. The most striking instances of gravitational lensing 
(e.g., multiple images, rings, arcs) are examples of rare phenomena caused by strong gravitational lenses, which greatly 
distort the images of distant galaxies. In contrast to this, weak gravitational lenses distort the images of distant galaxies 
very little but produce a net coherent pattern of image distortions in which there is a slight preference for the lensed 
galaxies to be oriented tangentially with respect to the direction vector that connects their centroids with the center 
of the gravitational potential of the lens. While weak lenses do not give rise to stunning individual images, they are 
detectable in a statistical sense via ensemble averages over many mildly-distorted images (e.g., |9], 1 10], 111], 1 12]). 

Provided the distance traveled by the light ray is very much greater than the scale size of the lens, it is valid to 
adopt the "thin lens approximation" in order to describe a gravitational lens. Consider a lens with an arbitrary 3- 
dimensional potential, <t>. In the thin lens approximation a conveniently scaled 2-dimensional potential for the lens 
(i.e., the 3-dimensional potential of the lens integrated along the optic axis) is given by 

= wk\l ®(P d e,z)dz, (2) 

D[D S c L J 

where 9 is the location of the lensed image on the sky, measured with respect to the optic axis, and D/ s , D/, and D s are 
angular diameter distances between the lens and source, observer and lens, and observer and source, respectively (e.g., 
(2l). It is then straightforward to relate the gravitational potential of the lens to the two fundamental quantities that 
characterize the lens: the convergence (k) and the shear (f). The convergence, which describes the isotropic focusing 
of light rays, is given by 

, , 1 /d 2 w d 2 w\ 



The shear describes tidal gravitational forces acting across a bundle of light rays and, therefore, the shear has both a 
magnitude, y = \ y\ + yf, and an orientation, (p. In terms of y/, the components of the shear are given by 



2\ddf del 

and 



= \ (^Z-^Z ) = 7(0)cos[2<p(0)] (4) 



(5) 



The effect of convergence and shear acting together in a gravitational lens is to distort the images of distant objects. 
Consider a source galaxy which is spherical in shape. In the absence of a gravitational lens, an observer would see an 
image of the galaxy which is truly circular. If a gravitational lens is interposed along the line of sight to the distant 
galaxy, the observer will see an image which, to first order, is elliptical and the major axis of the ellipse will be oriented 
tangentially with respect to the direction vector on the sky that connects the centroids of the image and the lens. That 
is, the circular source is distorted into an ellipse, and to first order the distortion consists of both a tangential stretch of 
(1 — K— y) -1 and a radial compression of (1 — fc+ y) _1 (e.g., (H). In the weak lensing regime, both the convergence 
and shear are small (k << 1 and y << 1). 

The fundamental premise in all attempts to detect weak lensing is that, in the absence of lensing, galaxy images have 
an intrinsically random ellipticity distribution. Gravitational lensing then introduces a shift in the ellipticity distribution 
that, in the mean, manifests as a tangential alignment of background sources around foreground lenses. The image of 
a distant galaxy can be approximated an ellipse with complex image ellipticity given by 

a 2 -b 2 

where a and b are the major and minor axes, respectively, and <p is the position angle. The complex image ellipticity is 
often referred to as the "image polarization" (e.g., 13911 ') and is computed in terms of flux-weighted second moments, 

Qij = ZhjWijXiXj, (7) 

hi 

where /,• ; is the intensity at a given pixel and Wij is a weighting function. The real and imaginary components of the 
image polarization are then given by: 

= 61.1-62,2 = 2Q h2 

' 61.1+62,2' 2 61,1+62.2' 

The observed image polarization for any one source is, of course, a combination of its intrinsic ellipticity and any 
ellipticity that is induced by lensing. In the limit of weak lensing, the observed image polarization, £ obs , is related 
to the intrinsic image polarization, e mt through a shift in the complex plane. Although we cannot determine £ lnt for 
any one particular source galaxy, we have that the mean intrinsic ellipticity distribution for an ensemble of source 
galaxies is (£ mt ) = since the galaxies should be randomly-oriented in the absence of lensing. An estimator for the 
shear induced by weak lensing is then y = (e obs ) /2 (e.g., |39]). This simple estimator does not reflect the fact that 
the way in which the shear alters the shape of a source depends u pon its intrinsic ellipticity, and in practice this is 
generally taken into account when computing the shear. See, e.g., 1 40], |41], and 1 42] for discussions of the "shear 
polarizability" and "shear responsivity" of sources. In addition, it is worth noting that, while it is common practice to 
approximate image shapes as ellipses, there will be some images that have been sufficiently distorted by galaxy-galaxy 
lensing that a mild bending, or "flexion", of the images will occur and such images cannot be accurately represented 
as ellipses. In principle, flexion of images can be used to detect weak lensing with a signal-to-noise that is increased 
over the common practice of fitting equivalent image ellipses |43], |44]. A preliminary application of this technique 
[43] has been carried out with the Deep Lens Survey [45], and it will be interesting to see how the technique is further 
developed and implemented in practice. 

The first attempts to detect systematic weak lensing of background galaxies by foreground galaxies (|46[, I47I0 were 
met with a certain degree of skepticism because the apparent distortion of the source galaxy images was rather smaller 
than one would expect based upon the typical rotation velocities of the disks of large spiral galaxies. The situation 





- 






1 1 1 1 1 | 1 Mill 1 Mill- 


2xl0- 3 








(b)-j 


lO- 3 










£ o 


Z I 


1 




- -1- t — i -$- m m • • • • • : 

( I $ ; 


-1x10-3 










-2x10-3 


- 






i i- 



0.1 1 10 100 

radius [arcrnin] 

FIGURE 2. a) Mean tangential shear computed about the lens centers in ~ 42 sq. deg. of the RCS 1 48]. Here foreground galaxies 
and background galaxies have been separated on the basis of apparent magnitude alone. Bright, lens galaxies have 19.5 < Rq < 21 
and faint, source galaxies have 21.5 < Rc < 24. b) Same as in a) except that here each background galaxy image has been rotated 
by 45°. This is a control statistic and in the absence of systematic errors it should be consistent with zero on all scales. Figure kindly 
provided by Henk Hoekstra. 

changed when Brainerd, Blandford & Smail 1 13] measured the orientations of 506 faint galaxies (23 < < 24) with 
respect to the locations of 439 bright galaxies (20 < r& < 23) and found that the orientation of the faint galaxies was 
inconsistent with a random distribution at the 99.9% confidence level. The faint galaxies showed a clear preference for 
tangential alignment with the direction vector on the sky that connected the centroids of the faint and bright galaxies, 
in agreement with the expectations of systematic weak lensing of the faint galaxies by the bright galaxies. 

Almost immediately, a number of similar investigations followed in the wake of Brainerd, Blandford & Smail lll3ll 
(149], |50], BUI, jH, |H, jH, 05]). These studies made use of a wide variety of data and analysis techniques, and 
all were broadly consistent with one another and with the results of Brainerd, Blandford & Smail II 311 (see, e.g., the 
review by Brainerd & Blandford lll4lo . The first truly undeniable detection of galaxy-galaxy lensing was obtained by 
Fischer et al. 15311 with 225 sq. deg. of early commissioning data from the SDSS, and it was this result in particular 
that helped to make the study of galaxy-galaxy lensing into a respectable endeavor, whereas previously many had 
considered the whole field rather dodgy at best. Fisher et al. l53ll demonstrated conclusively that even in the limit of 
somewhat poor imaging quality, including the presence of an anisotropic point spread function due to drift scanning, 
galaxy-galaxy lensing can be detected with very high significance in wide-field imaging surveys. In the last few years, 
detections of galaxy-galaxy lensing and the use of the signal to constrain the dark matter halos of field galaxies has 
improved dramatically (0], El, El, JH, JH, |H1, 0l, jH, IH) owing to a number of factors that include such 
things as very large survey areas, sophisticated methods for correcting image shapes due to anisotropic and spatially- 
varying point spread functions, and the use of distance information for large numbers foreground lens galaxies in the 
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FIGURE 3. Mean excess projected mass density around weak galaxy lenses in the SDSS 1 42]. Here ~ 1.27 x 10 5 lenses with 
spectroscopic redshifts and ~ 9.0 x 10 9 sources with photometric redshifts have been used in the calculation. The values of AL(r p ) 
shown in this figure have been corrected for the clustering of the sources around the lenses. Data kindly provided by Erin Sheldon. 



form of either spectroscopic or photometric redshifts. 

Figure |2] shows one example of the high statistical significance with which weak lensing due to galaxies is now 
being routinely detected. The result comes from an analysis of the distortion of the images of - 1.5 x 10 6 source 
galaxies due to — 1 .2 x 10 5 lens galaxies in the RCS 14811 . where the lens and source populations were separated solely 
on the basis of their apparent magnitudes. The top panel of Figure|2]shows the mean tangential shear computed about 
the lens centers which, because of the clustering of the lens galaxies, is not simply interpreted as the tangential shear 
due to individual lens centers. Instead, it is a projected (i.e., 2-dimensional) galaxy-mass cross-correlation function, 
and in order to compute the average properties of the halos of the lens galaxies it is necessary to, e.g., make use 
of Monte Carlo simulations that include all of the multiple weak deflections that the sources have undergone. The 
bottom panel of Figure |2 shows a control statistic in which the tangential shear about the lens centers is computed 
after rotating the images of the sources by 45°. If the signal in the top panel of Figure |2]is caused by gravitational 
lensing, the control statistic in the bottom panel of Figure[2]should be consistent with zero (and indeed it is). Note that, 
although the tangential shear about the RCS lenses persists to scales of order 0.5°, the shear on such large scales is 
not indicative of the masses of individual lens galaxies; rather it reflects the intrinsic clustering of the lenses. It is also 
worth noting that less than decade ago observers were struggling to measure a tangential shear of < 0.01 with a modest 
degree of confidence. Now, however, confident detection of tangential shears of < 0.0001 is effectively "routine" in 
these extremely large data sets. 

The mean tangential shear, fr(r p ), in an annulus of projected radius r p is related to the projected surface mass 
density of the lens through 

L c yr(r p ) = t(r<r p )-L(r p )=AL(r p ), (9) 

where E(r < r p ) is the mean surface mass density interior to the projected radius r p , T.(r p ) is the projected surface 
mass density at radius r p (e.g., Ir32ll . Ir33ll . O), and E c is the so-called critical surface mass density: 

E ^^k< (10) 



where c is the velocity of light and D s , D\, and D[ s are again angular diameter distances 0], The quantity A£(r p ) 
above is, therefore, a mean excess projected mass density. Shown in Figure [5] is the mean excess projected surface 
mass density in physical units of h Mq pc~ 2 for ~ 1.27 x 10 5 lens galaxies in the SDSS for which spectroscopic 
redshifts are known 1 42] . In addition to spectroscopic redshifts for the lenses, photometric redshifts were used for 
~ 9.0 x 10 9 source galaxies. Moreover, because the redshifts of the lens galaxies are known, A£(r p ) can be computed 
as a function of the physical projected radius at the redshift of the lens (rather than an angular scale). In Figure [3] 
AE(r p ) has been corrected for the clustering of the sources around the lenses via a function which is effectively a 
weighted cross-correlation function between the lenses and sources |42|. 

Having obtained a measurement of yr(8), or equivalently A£(r p ), constraints can then be placed on the nature of 
the dark matter halos of the lens galaxies by modeling the observed signal. As mentioned earlier, quite a bit of care 
has to be taken in doing this if the goal is to constrain the halo parameters as a function of, say, the host luminosity, 
color, or morphology (see, e.g., 11611 ). In the past few years, however, good constraints on the mass of an "average" 
halo associated with an L* galaxy, as well as fundamental differences between the halos of L* ellipticals versus L* 
spirals, have emerged from galaxy-galaxy lensing studies and it is those studies which are summarized below. 



PROBING HALO POTENTIALS WITH SATELLITE GALAXIES 

In order to use satellite galaxies to probe the potentials of host galaxies, one needs to define an appropriate sample 
of host and satellite galaxies. Unlike cosmology simulators who are blessed with full 6-dimensional phase space 
information, observers are, of course, limited to 3 dimensions (RA, DEC, and redshift). Given this limited information, 
then, one must base the selection criteria on projected radii (evaluated at the redshift of the host) and relative radial 
velocities, dv, of the candidate hosts and satellites. To guarantee that the dynamics of the satellites are determined 
solely by their host galaxy, the hosts must be determined to be "isolated" in some sense. That is, if another large, 
bright galaxy is too close to a candidate host galaxy to guarantee that the satellite orbits are affected solely by the 
candidate host, that candidate host and its satellites are rejected from the sample. Satellites must, necessarily, be 
fainter than their host, be found within some reasonable projected radius of the host, and have some reasonable line of 
sight velocity with respect to the host. 

There are a number of different selection criteria that have been used in the recent literature, and three sets of 
selection criteria that have been used in more than one investigation are summarized below: 

1 . Hosts must be at least 8 times brighter than any other galaxy that is within r p < 500 kpc and \dv\ < 1000 km sec . 
In addition, hosts must be at least 2 times brighter than any other galaxy that is within r p < 1 Mpc and 
\dv\ < 1000 km sec . Satellites must be at least 8 times fainter than their host, must be found within r p < 500 kpc, 
and must have \dv\ < 500 km sec -1 . Here h = 0.7 has been adopted (|20|, l64ll ). 

2. Hosts must be at least 2 times brighter than any other galaxy that falls within r p < 2.86 Mpc and \dv\ < 
1000 km s -1 . Satellites must be at least 4 times fainter than their host, must be found within r p < 714 kpc 
and must have \dv\ < 1000 km s -1 . Here h = 0.7 has been adopted (jH, (H, S, jH). 

3. Hosts must be at least 2.5 times brighter than any other galaxy that is within a projected radius of r p < 700 kpc and 
a relative radial velocity difference of \dv\ < 1000 km sec - 1 . Satellites must be at least 6.25 times fainter than their 
host, must be found within r p < 500 kpc, and the host-satellite velocity difference must be \dv\ < 500 km sec -1 . 
Here h — 0.7 has been adopted (|64], 1 68]). 

Although the above criteria may seem lax or even somewhat arbitrary, in the case of the first two sets of criteria, 
both the Milky Way and M31 would be excluded from the sample of hosts. That is, these particular selection criteria 
give rise to samples of unusually isolated host galaxies. In addition, both Prada et al. llzl and Brainerd l64ll adopted 
a number of different selection criteria in their investigations of the satellites of SDSS galaxies and concluded that 
there were no statistical differences between results that were obtained with different selection criteria. In other words, 
provided sufficiently "reasonable" criteria are adopted for selecting isolated hosts and their satellites, the results of the 
investigations are stable to modest differences in the details of those selection criteria. 

No matter what selection criteria are adopted, however, there will always be "interlopers" in the satellite data. 
Interlopers are galaxies that are falsely identified as satellites; that is, they pass the formal selection criteria, but they 
are not, in fact, dynamically associated with the host galaxy. The presence of interlopers will artificially inflate any 
measurement of the velocity dispersion of genuine satellites, and recent investigations of satellite dynamics (|27], 
1 65], 1 66], 1 67]) have corrected for the effects of interlopers by modeling the distribution of host-satellite velocity 
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FIGURE 4. Points with error bars show the observed distribution of velocity differences, ,/V(|dv|), for a subset of host-satellite 
systems in the 2dFGRS for which the host morphologies have been visually classified. Solid lines show the best-fitting "Gaussian 
plus offset" function, from which the velocity dispersion of the satellites, cv, and the fraction of interlopers, f;, is determined. Left 
panels: late-type hosts. Right panels: early-type hosts. Top panels: satellites located close to the host in projection on the sky. 
Bottom panels: satellites located far from the host in projection on the sky. A substantially larger value of o"V is obtained for the 
satellites of early-type hosts than for the satellites of late-type hosts. Note, too, that the fraction of interlopers increases significantly 
with the projected radius, r, of the satellites. 



differences as the sum of a Gaussian distribution (due to the genuine satellites) and a constant offset (due to the 
interlopers). Prada et al. l27ll used numerical simulations to show that this is a sensible way in which to correct for the 
effects of interlopers. Moreover, both Brainerd & Specian 1 66] and Prada et al. I27II have pointed out that an accurate 
determination of the velocity dispersion profile, o v (r p ), for satellite galaxies depends on a proper determination of the 
interloper fraction as an explicit function of the projected radius. That is, by purely geometrical effects, the interloper 
fraction is necessarily an increasing function of r p . An example of fitting a "Gaussian plus offset" to the distribution 
of velocity differences for late-type galaxies and early-type galaxies in the 2dFGRS is shown in Figure |4] One can 
clearly see from this figure that the velocity dispersion of the satellites is a function of the morphology of the host 
galaxy (being larger for early-type hosts than late-type hosts), and that the interloper fraction increases with projected 
radius. 

The above "Gaussian plus offset" fit to the distribution of host-satellite velocity differences accounts for the fact 



that the number of interlopers is a function of projected radius, and it assumes a priori that the number of interlopers at 
a given projected radius is constant with \dv\. Recently, however, van den Bosch et al. llfwll used simulations of galaxy 
redshift surveys to investigate this and found a sharp increase in the number of interlopers for small relative velocities, 
van den Bosch et al. 1 69] note, however, that the value of a v that is determined from a simple "Gaussian plus offset" 
fit is not strongly affected by the fact that the number of interlopers varies with \dv\. This is because the best-fitting 
value of (7 V is rather insensitive to the precise value of the interloper fraction. Brainerd |67] also finds that the number 
of interlopers is larger for small values of |<jy| than it is for large values of \dv\, but that the effect is not nearly as 
pronounced as found by van den Bosch et al. 1 69] . Given the size of the error bars on the distribution of host-satellite 
velocity differences in the current observational samples, then, it would appear that the simple "Gaussian plus offset" 
fit to the distribution of velocity differences is more than adequate to the task of estimating o v (r p ). 



THEORY: "UNIVERSAL" (NFW) HALOS VS. ISOTHERMAL HALOS 

High-resolution CDM simulations have established the existence of a "universal" density profile for dark matter 
halos which results from generic dissipationless collapse (e.g., H, H, 0], 0, 0, jZJ, 0, 0, 0)- This 
density profile fits objects that span roughly 9 orders of magnitude in mass (ranging from the masses of globular star 
clusters to the masses of large galaxy clusters) and applies to physical scales that are less than the "virial" radius, r2oo- 
Conventionally, r2oo is defined to be the radius at which the spherically-averaged mass density reaches 200 times the 
critical mass density (e.g., |0L &). 

Navarro, Frenk & White |5], |6], |7] showed that the universal density profile for dark matter halos was fitted well 
by a function of the form 

P« = — 7, (ID 

(r/r s )(l+r/r s ) 2 

and halos having such a density profile are generally referred to as "NFW" halos. Here p c = - &7C £' is the critical 
density of the universe at the redshift, z, of the halo, H(z) is Hubble's parameter at that same redshift, and G is 
Newton's constant. The scale radius r s = r2oo /c is a characteristic radius at which the density profile agrees with the 
isothermal profile (i.e., p (r) r~ 2 ), c here is a dimensionless number known as the concentration parameter, and 

_ 200 c 3 
c ' = ~ln(l+c)-c/(l+c) (12) 

is a characteristic overdensity for the halo. 

Formally, the above fitting function for the radial density profiles of CDM halos converges to a steep, cuspy profile: 
p(r) r . The NFW fitting formula, however, was never intended to be extrapolated to very small radii (i.e., radii 
smaller than the practical resolution limits of the simulations) and much fuss has been made over whether observed 
galaxies actually show such cuspy inner density profiles (e.g., 0, 0, 0, 0, 0, 0, 0, 0). More 
recent numerical work has shown that the density profiles of CDM halos do not, in fact, converge to a well-defined 
asymptotic inner slope (e.g., 1 84], 1 85] 1 86], 1 87]), and it has become increasingly clear that fair and direct comparisons 
of simulated galaxies with observed galaxies on very small physical scales is an extremely challenging thing to do (e.g., 
0,0). 

Weak lensing and satellite dynamics do not have the ability to provide any information whatsoever on the cuspiness 
(or lack thereof) in the central regions of galaxies. Instead, both are governed by the large-scale properties of the halos 
(i.e., the regime in which the NFW profile is known to be an excellent description of the density profiles of CDM halos) 
and, at least in principle, both have the potential to discriminate between NFW halos and simpler singular isothermal 
sphere halos. 

The radial density profile of a singular isothermal sphere halo is given by 

P« = 2^ (13) 

(e.g., 1 90]), where <J V is the velocity dispersion. The isothermal sphere is characterized by the single parameter d v , 
which is constant as a function of radius. A key prediction for NFW halos, however, is that the radial velocity 
dispersion will have a strong dependence upon the radius and this, of course, is inconsistent with the constant value of 



the velocity dispersion that characterizes singular isothermal spheres. Specifically, on sufficiently small scales aJr) 
should increase with radius, and on large scales <J r (r) should decrease with radius. Hoeft, Miicket & Gottlober Il9l1l 
have shown that the radial velocity dispersion of NFW halos can be fitted by a function of the gravitational potential, 
<t>(r), of the form: 

K"l 1/2 
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Note that G r (r) above is not the "line of sight" velocity dispersion, since r is a true 3-dimensional radius in eqn. MAI . 
The parameters a and K have values of a — 0.29 ± 0.04 and K — 0.41 ± 0.03, and <t> out is the outer potential of the 
halo. Therefore, we expect the dynamics within an NFW halo to differ fundamentally from the dynamics within an 
isothermal sphere halo. 

In the case of weak lensing, NFW halos give rise to a distortion in the images of distant galaxies that differs 
somewhat from the distortion that would be yielded by an isothermal sphere halo (e.g., ll92ll . l93lo . The radial 
dependence of the shear for the isothermal sphere is given by: 
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(e.g., 1 81]). Here c is the velocity of light and D s , Dj, and D/ s are again angular diameter distances. In the case of NFW 
halos, the radial dependence of the shear is given by: 
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where x = r p /r s , £ c is the critical mass density for gravitational lensing given by eqn. (1101 . and the functions g<,> (x) 
are explicitly independent of the cosmology: 
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(e.g., |9J, 19J). 

In the following sections I summarize the most recent attempts to study the dark matter halos of field galaxies 
through satellite dynamics and weak lensing, including attempts to distinguish between isothermal and NFW potentials 
on the basis of the velocity dispersion profile and on the weak lensing shear. 



OBSERVED VELOCITY DISPERSION PROFILES 

At best, galaxy-galaxy lensing and satellite dynamics have the potential to constrain the dependence of the line of 
sight velocity dispersion on the projected radius, <7 v (r p ). Determining O v {r p ) has proven to be quite a challenge to 
galaxy-galaxy lensing studies, in large part because the shear profiles of NFW lenses and isothermal sphere lenses are 
not dramatically different, except on the very smallest (r < r s ) and very largest (r > r v ; r ) scales 1 93 ] . To date, only one 
tentative measurement of O v {r p ) has been made from observations of galaxy-galaxy lensing 1 24] . Kleinheinrich et al. 
1 24] modeled the lens galaxies in the COMBO-17 survey as singular isothermal spheres with velocity dispersions that 
scaled with luminosity as 
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FIGURE 5. Velocity dispersion profiles for satellites of SDSS host galaxies [27]. Circles: host galaxies with —20.5 < Mb < 
—21.5, Squares: host galaxies with —19.5 < Mb < —20.5. Left panel: "raw" velocity dispersion profiles prior to correction for 
contamination by interlopers. Right panel: velocity dispersion profiles after correction for contamination by interlopers. After 
correction for interlopers, O v (r p ) for the satellites of the fainter hosts is consistent with the expectations for an NFW halo with 
M200 = 1-5 x 10 12 M e , and O v (r p ) for the satellites of the brighter hosts is consistent with the expectations for an NFW halo with 
M200 = 6 x 1O I2 M . Here h = 0.7 has been adopted. 



where a* is the line of sight velocity dispersion of the halo of an L* galaxy. Kleinheinrich et al. 1I24I1 fixed T] to be 
0.35 and determined best-fitting values of a* for projected radii in the range 20 hr x kpc < r p < r max . When they 
considered all lenses in their sample, Kleinheinrich et al. 12411 found a* ~ 139 km sec _1 for r max = 50 h kpc, 
(7* ~ 1 64 km sec for r ma x = 150 h- 1 kpc, and a* ~ 123 km sec for r max = 500 h 1 kpc. This suggests a velocity 
dispersion profile that rises at small radii, reaches a maximum, then decreases at large radii. However, the formal error 
bars on these measurements show that all of these values of a* agree to within one to two standard deviations. In 
addition, it should be kept in mind that each of these measurements of a* is not independent (as they would be if 
a differential measurement of <T*(r p ) were made), so the data points and their error bars are all correlated with one 
another. 

Considerably stronger constraints on the dependence of the halo velocity dispersion with projected radius have come 
from the most recent investigations of the motions of satellites about host galaxies. In particular, both Prada et al. 12711 
and Brainerd 16711 have measured decreasing velocity dispersion profiles for the satellites of host galaxies in the SDSS 
and 2dFGRS, respectively. Although they used different data sets and different host-satellite selection criteria, both 
Prada et al. 1 27] and Brainerd 1 67] used the same technique to make measurements of the velocity dispersion profiles. 
That is, the distribution of velocity differences, N(\dv\), for satellites found within projected radii of r m j n < r p < r max 
was modeled as a combination of a Gaussian and an offset due to interlopers. In both studies, the interloper fraction 
was determined separately for each of the independent radial bins. 

Prior to correcting for the contamination of interlopers, Prada et al. 1I27I1 found a velocity dispersion profile, o v {r„), 
that increased with projected radius. After the removal of the interlopers, however, Prada et al. l27ll found decreasing 
velocity dispersion profiles in both cases. The corresponding velocity dispersion profiles are shown in Figure [5] 
Moreover, their corrected velocity dispersion profiles were fitted well by the velocity dispersion profiles of NFW 
halos with virial masses of 1.5 x 10 12 M Q (hosts with absolute magnitudes —19.5 < M B < —20.5) and 6 x 10 12 M Q 
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FIGURE 6. Velocity dispersion profiles for satellites in the final data release of the 2dFGRS and the flat, A-dominated GIF 
simulation 1 67]. Here h = 0.7 has been adopted. 



(hosts with absolute magnitudes —20.5 < Mg < —21.5). Since Prada et al. II27T adopted a value of h = 0.7 and since 
the absolute magnitude of an L* galaxy is M* B ~ —19.5, these results suggest that the virial mass of the halo of an L* 
galaxy is < 10 x IO^^Mq. 

Brainerd |67) selected hosts and satellites from the final data release of the 2dFGRS using criteria identical to 
those of Sample 3 in Prada et al. l27ll . In addition, she used these same criteria to select hosts and satellites from the 
present epoch galaxy catalogs of the flat, A-dominated the GIF simulation 1 94] . This is a publicly-available simulation 
which includes semi-analytic galaxy formation in a CDM universe. Brainerd 1 67] restricted her analysis to hosts with 
luminosities in the range 0.5 L* < L < 5.5 L*, and found a roughly similar number of hosts and satellites in both the 
2dFGRS (1345 hosts, 2475 satellites) and the GIF simulation (~ 1200 hosts, ~ 4100 satellites, depending upon the 
viewing angle). Like Prada et al. l27ll . Brainerd 1 67] obtained a decreasing velocity dispersion profile for the satellites 
of the 2dFGRS galaxies once the effects of interlopers were removed. In addition, excellent agreement between o v (r p ) 
for the 2dFGRS galaxies and a v (r p ) for the GIF galaxies was found, showing consistency between the motions of 
satellites in the 2dFGRS and the expectations of a A-dominated CDM universe. See Figure|6] 

Further, Brainerd [67 ] divided her sample of 2dFGRS host galaxies into thirds based upon the spectral index pa- 
rameter, 77 1 95], and computed the dependence of the velocity dispersion profile on host spectral type. The subsamples 
corresponded to hosts which are expected to have morphologies that are approximately: (i) E/S0, (ii) Sa, and (iii) 
Sb/Scd. The median luminosities of the hosts in the subsamples were all fairly similar: (i) 2.64 Lt , (ii) 2.25 L* b , and 
(iii) 2.11 L* bj . The velocity dispersion profiles of all three samples decreased with radius and, moreover, O v {r p ) was 
found to have a much higher amplitude and steeper decline for the satellites of early-type hosts than it did for the 
satellites of late-type hosts. See Figure0 Although there is some difference in the median luminosities of the hosts in 
the subsamples, the difference is too small to have a significant effect on the velocity dispersion profiles. Therefore, the 
results of Brainerd 1 67] seem to indicate that early-type galaxies have deeper potential wells (and hence more massive 
halos) than late-type galaxies. 

Previous work on the dependence of a, with projected radius using SDSS galaxies i65Tl and 2dFGRS galaxies l67il 
concluded that o v (r p ) was consistent with an isothermal profile; i.e., O v {r p ) = constant. In both of these investigations, 
the hosts and satellites were selected in a manner that was identical to that of Sample 3 in Prada et al. 12711 . In both 
previous analyses, however, the number of hosts and satellites was significantly smaller than the more recent studies, 
and the formal error bars were correspondingly larger. In addition, the original analysis of SDSS host-satellite systems 
1 65] neglected to account for the fact that the interloper fraction increases with radius, which would have biased 
measurements of a v at large r p towards values which are higher than the actual satellite velocity dispersion at those 
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FIGURE 7. Velocity dispersion profiles for satellites in the final data release of the 2dFGRS as a function of the host spectral 
parameter, r\ 1 67]. The morphology of the hosts is expected to be roughly E/SO in the left panel, Sa in the middle panel, and Sb/Scd 
in the right panel. The median luminosities of the subsamples in each of the panels is somewhat different, but the difference is too 
small to account for the differences in the velocity dispersion profiles. Here h = 0.7 has been adopted. 



radii. 

Even more recently, Conroy et al. ll%ll used satellites of z ~ 0.8 host galaxies in the DEEP2 survey to investigate 
O v (r p ). DEEP2 (Deep Extragalactic Evolutionary Probe 2) is being carried out with the DEIMOS spectrograph at 
the Keck-II telescope, and will ultimately collect spectra of ~ 60,000 galaxies with redshifts of 0.7 < z < 1.4 to 
a limiting magnitude of R^g = 24.1 19 711 . Unfortunately, the survey is still far from complete and only 61 isolated 
host galaxies (having a total of 75 satellites) were found in the current DEEP2 data. Because of this, the errors on 
o v (r p ) are large, and formally <7 v (r p ) for the DEEP2 galaxies is fitted well by a constant value: <7 V (110 h kpc) = 
162+44 km sec" 1 , a v (230 ft -1 kpc) = 136+|^ km sec" 1 , a v (320 hr 1 kpc) = 155+|| km sec" 1 . Therefore, isothermal 
halos for the DEEP2 galaxies cannot be ruled out at the moment. Conroy et al. ll96ll show, however, that their 
velocity dispersion measurements are consistent with expectations for NFW halos with virial masses in the range 
3.5 x 10 12 hr^ Mq < M200 < 8.0 x 10 12 h~ x M . This is in good general agreement with the results of Prada et al. 
ll27ll . especially considering that the DEEP2 hosts are of order one magnitude brighter than the SDSS hosts (i.e., the 
virial mass implied for the halos of the brightest galaxies in the SDSS sample is ~ 4 x 10 l2 /i _1 M Q ). At the moment, 
however, the DEEP2 data are too severely limited by small number statistics to place strong constraints on the nature 
of the dark matter halos of galaxies with redshifts of order unity. 



HALO MASSES AND GALAXY MASS-TO-LIGHT RATIOS 

Although it ought to be straightforward and even easy to compare the halo masses and galaxy mass-to-light ratios 
that are obtained from different studies, in practice it is rather like comparing persimmons to tomatoes; i.e., they are 
vaguely similar on the inside and outside, but they are definitely not interchangeable. The fundamental problem is that 
it is simply not possible to measure the "total" mass of a galaxy halo (since it is not possible to say where such a halo 
"ends") and, hence, all halo masses are simply masses that are contained within some physical radius of the center 
of the halo. Along those same lines, and given that velocity dispersion profiles of NFW halos decrease with radius, 
if one wants to compare the results of two investigations which have measured a velocity dispersion averaged over 
some large scale, it is important that those scales be identical. That is, suppose a single measurement of <7 V is made 
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FIGURE 8. Isothermal sphere models for the galaxy-galaxy data from COMBO-17 1 60]. Joint constraints (Id, 2cx, and 3a) on 
the velocity dispersion, <T* , of the halos of L* galaxies and the index of the Tully-Fisher/Faber-Jackson relation, rj. Here the weak 
lensing signal has been averaged over scales r p < 150 h~ l kpc. Left panels: all lenses, cr* = 156+ |g km sec -1 , tj = 0.28+QQg. Right 
panels: red lenses (2579 galaxies, cr* red = I8O+3Q km sec -1 ) Left panels: blue lenses (9898 galaxies, <7* blue = 126+^ km sec -1 ). 
Figure kindly provided by Martina Kleinheinrich. 

by averaging over scales r < 100 /i -1 kpc in one study and a single measurement of a v is made by averaging over 
scales r < 200 hT x kpc in another. If the second measurement of <7 V is lower than the first by some significant amount, 
that does not necessarily mean that the values are in disagreement. They would be in disagreement if both halos were 
isothermal spheres, but if the halos are NFW objects, then it is only to be expected that the second measurement would 
be lower than the first. 

A more subtle problem is the definition of the "virial radius" in the context of NFW halos. While r2oo was originally 
proposed as the radius at which the interior mass density is 200 times the critical mass density (e.g., Q], 0,10]), it 
is not at all uncommon to find that investigators who have fit NFW models to their data have defined the virial radius 
as the radius at which the interior mass density is 200 times the mean mass density of the universe. Therefore, what is 
meant by a "virial mass" in the context of an NFW fit to data can (and does) vary from investigation to investigation, 
and a certain amount of care has to be taken when comparing such results. Despite the difficulties of comparing the 
conclusions of different studies, I will forge ahead because it is becoming clear that a consistent picture really is 
emerging on the topic of the masses of the halos of field galaxies, and their corresponding mass-to-light ratios. The 
weak lensing studies yield results that are by and large consistent with each other, and the dynamical studies seem 
to be in general agreement with the trends in the weak lensing data: the halos have masses that are consistent with 
expectations for galaxy-sized halos in CDM, and there are real, physical differences between halos surrounding (i) 
early-type and late-type galaxies and (ii) high-luminosity and low-luminosity galaxies. 

M and M/L from Galaxy-Galaxy Lensing 

In the case of galaxy-galaxy lensing, it is not possible at the moment to discriminate between shear profiles that 
are caused by NFW versus isothermal galaxy halos. Therefore, investigators will often choose one or the other to 
constrain the properties of the halos that are producing the lensing signal. In the case of isothermal sphere halos, the 
velocity dispersions of the lens galaxies used to model the observed signal are often chosen to scale as in eqn. \\9\ 
above, (ov/cr*) = (L/L*)* 1 , where again a v is the velocity dispersion of a halo that contains a galaxy of luminosity 
L, and <7* is the velocity dispersion of the halo of an L* galaxy. Hoekstra et al. |48] used this approach with their 
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FIGURE 9. Constraints on the circular velocity at r = r20o and the scale radius, r s , for lenses in the RCS that have been 
modeled as having NFW-type halos Formally, the best-fitting values of the circular velocity, scale radius and virial mass are: 
^200 = 162 ±8 km sec -1 , r s = 16.2+| \ hr l kpc, and M200 = 8.4 ± 1.1 x 10 11 h~ l M . Here r2oo is defined as the radius at which 
the mean interior mass density of the halo is equal to 200p c . The dashed line shows the predictions of the NFW theory, in which 
V200 and ^200 are not independent parameters. Figure kindly provided by Henk Hoekstra. 



RCS data, as did Kleinheinrich et al. 1 60] with their COMBO-17 data. When all lenses and sources were used 
in the investigations, and when the lensing signal was averaged over an identical scale (r < 350 h~ l kpc), both 
the RCS and COMBO-17 results are in very good agreement with each other. In particular, Hoekstra et al. |48] 
find c* = 136 ± 8 km sec for an adopted value of 77 = 0.3, and Kleinheinrich et al. 0] find a* = 138+i| and 
77 = 0.34+ 

Em- Further, Kleinheinrich et al. [60] find that there are clear differences in the halos surrounding "blue" 
galaxies (rest frame colors of (U — V) < 1.15 — 0.3 lz — 0.08 [My — 5 log/i + 20]) and those surrounding "red" galaxies 
(the remainder of the sample). That is, the red COMBO-17 lens galaxies have a higher velocity dispersion than the 
blue COMBO-17 lens galaxies, but both have a similar value of the index 77 above. See Figure|H] 

In addition, Guzik & Seljak 1 16], Hoekstra et al. |48], and Kleinheinrich et al. 1 60] have all used NFW halos to model 
their lens galaxies, and all find very reasonable fits to their lensing signals. Further, the derived values of the NFW 
virial masses of the halos of L* galaxies are in quite good agreement amongst these studies when they are determined 
in similar band passes (e.g., r) and with identical definitions of the virial radius 1 60 ] : M* ir = 8.96 ± 1.59 x 10 /z 1 M Q 
^,M* k = 8.4±0.7 x lO 11 /!" 1 M i48ll.andM* ir = 7.8 + H x 10 11 /i" 1 M [60]. These are also in remarkably good 
agreement with the virial mass implied for the halos of L* galaxies by the dynamical analysis of Prada et al. 12711 
(e.g., M* iT - 10 x I0 n h- 1 M & ). Shown in Figure are 1(7, 2(7 and 3(7 confidence limits on a joint-parameter fit of 
the circular velocity at r2oo, V200, and scale radius, r s , for the lenses in the RCS data ll48ll . Note that in the analysis 
of the RCS data, V200 and r s were allowed to vary freely, while, to within some scatter, these parameters are strongly 
correlated in the NFW theory (i.e., the NFW model is in essence specified by a single parameter). The dashed line 
in Figure |9] therefore shows the prediction for a strict adherence to the NFW theory (i.e., V200 and r„ are correlated 
appropriately), and the fact that the theoretical NFW line passes so well through the contours gives a certain amount of 
confidence that the NFW model is a very good fit to the data. Kleinheinrich et al. 1 60] find good fits of the NFW model 
to their data and, moreover, find that both the virial radii of the halos and the parameter tj are dependent upon the rest 
frame colors of the galaxies, with red galaxies having a somewhat larger virial radius (and, hence, larger virial mass) 
than blue galaxies. See Figure[K)] Here tj is defined not as in eqn. dl9> , since the velocity dispersion is a function of 
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FIGURE 10. NFW halo models of the galaxy-galaxy lensing data from COMBO-17 |H3l- Joint constraints (lcr, 2a, and 3a) 



on tj and the virial radii of the halos of L* galaxies are shown. Left panel: all lenses, tj = 0.30 
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projected radius in the NFW model, but rather it is defined as: 
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in analogy to the Tully-Fisher and Faber-Jackson relations (see 1 60]). In this case, r* ir is the virial radius of the halo 
of an L* galaxy, defined at 200 times the mean mass density of the universe. The variation of tj with galaxy color and 
its implications for the mass-to-light ratios of the galaxies will be discussed below. 

A particularly detailed study of the masses of lensing galaxies as a function of their color was carried out by Guzik 
& Seljak [ 16] for - 3.5 x 10 4 lenses and - 3.6 x 10 6 sources in the SDSS. All of the lens galaxies have spectroscopic 
redshifts in this case, and all of the halos were modeled as NFW objects in the context of the "halo model". In all 5 
of the SDSS band passes, Guzik & Seljak 1 16] find that the virial masses of L* ellipticals exceed those of L* spirals 
though, unsurprisingly, the amount by which the masses of the ellipticals exceeds those of the spirals is a strong 
function of the band pass. In the redder bands, the masses of the ellipticals exceed those of the spirals by a factor of 
~ 2 to ~ 2.5, while in g' the difference is a factor of ~ 6 and in u' the difference is close to an order of magnitude. 
Although it is difficult to make direct comparisons between the two studies (because of the differing definitions of 
the virial radius and the different definitions of the subsamples of galaxies), there is good general agreement between 
the results of Guzik & Seljak 1 16] and Kleinheinrich et al. 16011 : when the galaxy-galaxy lensing signal is detected 
red band passes (e.g., R, r 1 ) and the lenses are modeled as NFW objects, the virial masses of red/early-type galaxies 
exceed those of blue/late-type galaxies by a factor of order 2. 

In addition to the halos of early-type lenses having more mass than those of late-type lenses, the weak lensing 
work of Sheldon et al. i42ll indicates that, again, in all 5 SDSS band passes, the projected excess surface mass density 
increases with the luminosity of the lens. Sheldon et al. I42ll separated their ~ 1.27 x 10 5 lenses into 3 magnitude 
bins (high, middle, and low luminosity), and the magnitude cuts differ for the different band passes. (See Table 2 
of Sheldon et al. 1 42] for a complete list of the magnitude cuts as a function of band pass.) In the case of the r' 
data, the "high" luminosity galaxies have a mean absolute magnitude of -22.5, the "middle" luminosity galaxies have 
a mean absolute magnitude of -21.9, and the "low" luminosity galaxies have a mean absolute magnitude of -20.5. 
These mean luminosities correspond roughly to 4.5L* ("high"), 2.7L* ("middle") and 0.8L* ("low") in the r* band. 



In all cases, AL(r„) for the "high" luminosity galaxies exceeds that of the "medium" and "low" luminosity galaxies, 

and for r p < 1 h Mpc, the difference corresponds to an approximately constant multiplicative factor. Specifically 

at r p ~ 100 hr l Mpc, however, AE for the high luminosity lenses in Sheldon et al. 1 42] exceeds that for the low 

luminosity lenses by a factors of ~ 3 in ~ 5 in g' , ~ 5 in r', ~ 7 in i', and ~ 7 in z! (e.g., Figure 14 of Sheldon et 

al. I42in . Similar trends (i.e., higher projected excess surface mass density for more luminous lenses) were found by 

Seljak et al. l6lll in their galaxy-galaxy lensing analysis of SDSS data. 

Lastly, although there is reasonable agreement regarding the relative increase in mass for the halos of early-type lens 

galaxies versus late-type lens galaxies at fixed luminosity (i.e., L*), there is some disagreement over the dependence 

of the mass-to-light ratio on the luminosity of the host. Specifically, in their redder bands Guzik & Seljak 11611 find 

that the mass-to-light ratio goes as M/L« ^o.4±o.2 f or ^ > j*^ SU gg es tive of a mass-to-light ratio that increases 

with luminosity. Kleinheinrich et al. 16011 . however, find that M/L for their sample of lenses is more consistent with 

-o io +0,48 m m 

a constant value: M/L °= L ' -0.36. Both Guzik & Seljak 1 16] and Kleinheinrich et al. 1 60] agree, however, that the 

mass-to-light ratio of red/early-type L* lens galaxies exceeds that of blue/late-type L* lens galaxies by a factor of 

~ 2 to ~ 2.5 in the redder bands. 



M and M/L from Satellite Dynamics 

In the 1 990's, Zaritsky et al. j20ll and Zaritsky & White 1 98 ] used the velocity differences between a small number of 
isolated spiral galaxies and their satellites to show that the halos of the spirals were massive and extended to large radii: 
M(150 h kpc) ~ 1 to 2 x 10 l2 h~ l M & . Moreover, Zaritsky et al. 12011 found a somewhat curious result: the velocity 
difference between their 115 satellites and 69 hosts was independent of the inclination corrected H-I line width of the 
host and was, therefore, independent of the luminosity of the host (through, e.g., the Tully-Fisher relation). At fixed 
large radius, then, this would imply that M/L for the spiral hosts decreased as M/L« L _1 . 

More recent investigations of halo masses and corresponding mass-to-light ratios from satellite dynamics have led 
to rather a large assortment of conclusions. McKay et al. |65] and Brainerd & Specian |66] used the dynamics of the 
satellites of SDSS galaxies and 2dFGRS galaxies, respectively, to constrain the dynamical masses of the halos of the 
host galaxies interior to a radius of r = 260 hr x kpc. Both used an isothermal mass estimator of the form 

<o = ^, (2D 

where a v is the line-of-sight velocity dispersion. Both felt this assumption was justified because both found that 
their velocity dispersion profiles were consistent with a constant value. In the case of McKay et al. 1 65], however, 
no correction for an increasing number of interlopers with projected radius was made and this may have led to an 
incorrect conclusion that O v (r p ) was independent of r p . In the case of Brainerd & Specian [66], the increasing number 
of interlopers at large r p was taken into account, but only galaxies from the 100k data release of the 2dFGRS were 
used (i.e., roughly half as many galaxies as in the full data release), and although <7 v (r p ) was consistent with a constant 
value in their data, the later analysis by Brainerd [67] showed that this was simply due to the rather large error bars 
in Brainerd & Specian 1 66] . This being the case, the mass-to-light ratios published by these two studies are suspect 
at some level, but it is unclear at the moment just how suspect they may actually be. That is, while it is true that the 
velocity dispersion profile of NFW halos decreases with radius, the fall-off in O v (r p ) is not particularly sharp and it is 
not obvious how badly isothermal mass estimates of the form in eqn. (12 It . which are based on an average value of a v , 
will compare to proper NFW mass estimates. 

Formally, McKay et al. ||65| found that in all 5 SDSS band passes, M^ /L was roughly constant for L > L* , and that 
the value of M^l/L was a strong function of the band pass (being systematically higher in the blue bands than in the red 

bands). Brainerd & Specian 1 66] found that for L > 2L*, M^Iq/L was a constant for dynamical analyses that included 
(i) all 809 hosts in their sample and (ii) 159 hosts that had been visually classified as early-type (E/S0). However, 
much like the results of Zartisky et al. 1 20], Brainerd & Specian 1 66 ] found that M^/L decreased as M^Iq/L °= L~' 
for 243 hosts that had been visually classified as spirals. This latter result remains puzzling, and is certainly in need of 
further investigation with larger data sets. 

In their analysis of the dynamics of the satellites of SDSS host galaxies, Prada et al. l27ll found that the velocity 
dispersion of the satellites scaled with host luminosity as ff,, °= L° 3 (i.e., in good agreement with the local fi-band 
Tully-Fisher relationship 1 99]) for satellites with projected radii r p < 120 kpc. (Recall, too, that in this study O v {r p ) 
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FIGURE 11. Dependence of satellite velocity dispersion on host absolute magnitude for SDSS galaxies 1 27]. Filled circles: a v 
computed using satellites with 20 kpc < r p < 120 kpc. Open circles: o",, computed using satellites with 250 kpc < r p < 350 kpc. For 
small projected radii the velocity dispersion scales as c>,> °= L 0,3 , in good agreement with the local B-band Tully-Fisher relationship. 
For large projected radii a v °= L 5 . Here h = 0.7 has been adopted. 



was specifically corrected for the increase in interlopers at large r p .) In addition, Prada et al. |27] found that for 
satellites at large projected radius, 250 kpc < r p < 350 kpc, the velocity dispersion scaled with luminosity as rj v oc L 5 
(i.e., steeper than expected from the Tully-Fisher relation). See FigurefTTI 

Similar to Prada et al. 12711 . Brainerd 1 67] also computed the dependence of the small-scale velocity dispersion of 
satellites on host luminosity. See Figure[21 Like Prada et al. 1 27], Brainerd |67] corrected for the fact that the interloper 
fraction is an increasing function of projected radius and overall, she found excellent agreement between the velocity 
dispersions of satellites with projected radii r p < 120 kpc in the 2dFGRS and GIF simulations. The velocity dispersions 
of the 2dFGRS satellites were, however, seen to scale with host luminosity as a v « L^ 45±0 10 , which is only marginally 
consistent with the results of Prada et al. l27ll and the local Z?-band Tully-Fisher relationship. 

Prada et al. [27] have shown (e.g., their Figure 12) that the dependence of the line of sight velocity dispersion 
on the virial mass of NFW halos scales as o Y M° ir 38 for the case that rj,, is computed as an average over scales 
20 kpc < r p < 100 kpc, and that rj,, °c M^ ir 50 for the case that <7,. is computed at r p ~ 350 kpc. Combining this with 
their results for the dependence of o v on L at different scales leads to the conclusion that on scales r p < 120 kpc, 
M v ij./L °< L~ - 2 while on scales r p ~ 300 kpc, M v ; r /L is a constant. Similarly, if the halos of the 2dFGRS galaxies 

J— l a 2+0-3 

studied by Brainerd 1 67] are assumed to be NFW objects, the implication is that M v i r /L °= L ' -° 1 for the 2dFGRS 
hosts (again, computed on scales r p < 120 kpc). 

While it certainly cannot be said that there is a consensus from weak lensing and satellite dynamics as to the exact 
dependence of the galaxy mass-to-light ratio on L, it does seem to be the case that all of these studies point towards a 
dependence of M v ; r /L on L that is, at most, rather weak. That is, with the notable exception of the Brainerd & Specian 
[66] result for late-type galaxies, all of the recent determinations of M/L for L>L* find that, to within 2a, M/L is 
independent of L. In addition, when the weak lenses and host galaxies are each modeled as NFW objects, a fairly 
consistent value of the average virial mass of the halos of L* galaxies is found: ~ (8 — 10) x 10 11 h~ l M Q . Further, 
it seems to be clear that both weak lensing and satellite dynamics indicate that the masses of the halos of early-type 
galaxies are larger than that of late-type galaxies, and that at fixed luminosity the mass-to-light ratios of early-type 
galaxies are larger than those of late-type galaxies. 
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FIGURE 12. Dependence of satellite velocity dispersion on host luminosity for satellites with projected radii r p < 120 kpc in 



both the 2dFGRS and the flat, A-dominated GIF simulation 1671. Dotted line shows o",, « L 
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NON-SPHERICAL HALOS 



Although the simple isothermal sphere can reproduce the flatness of the rotation curves of the disks of spiral galaxies 
at large radii, there are both observational and theoretical arguments in favor of halos which are flattened, rather 
than spherical. Direct observational evidence for halo flattening that has come from studies of individual galaxies is 
somewhat scarce, however, owing to the fact that there are relatively few galaxies for which the shape of the halo 
potential can be probed directly. Nevertheless, the evidence for flattened halos of individual galaxies is diverse and 
includes such observations as the dynamics of polar ring galaxies, the geometry of X-ray isophotes, the flaring of 
HI gas in spirals, the evolution of gaseous warps, and the kinematics of Population II stars in our own Galaxy. In 
particular, studies of disk systems which probe distances of order 15 kpc from the galactic planes suggest that the 
ratio of shor test to longest principle axes of the halos is c/a = 0.5 ±0.2 (see, e.g., the comprehensive review by 
Sackett [ 100] and references therein). Studies of a number of strong lens galaxies hav e also suggested that the mass 
distributions of the lenses are not precisely spherical. For example, Mailer et al. 1 101] found that, provided the disk 
mass is small compared to the halo mass, the halo of the spiral galaxy which lenses the quasar B 1600+434 is consistent 
with c/a = 0.53. In addition, the 17 strong lens systems studied by Keeton, Kochanek & Falco 1 102] showed some 
preference for flattened mass distributions, although extremely flattened (i.e., "disky") mass distributions were ruled 
out. Finally, a recent analysis of the luminous halos of 1047 edge-on disk galaxies in the SDSS suggests that the old 
stellar populations of these galaxies consist of moderately flattened spheroids with axis ratios of c/a ~ 0.6 1 103]. 

On the theoretical side, high-resolution simulations of dissipationless CDM models consistently produce markedly 
non-spherical galaxy halos with a mean projected ellipticity of e ~ 0.3 (see, e.g., 1 104], 1 105|). It is known, however, 
that the dark matter will react to the condensation of baryons during galaxy formation (e.g., Il06l0 and that the res ultin g 
increase in the central density leads to a more spherical shape than if dissipation were not considered (e.g., Col). 
Recent simulations performed by Kazantzidis et al. 1 108] show that on scales r « r ya , the effects of gas cooling 
cause a substantial circularization of the mass density profile, leading to a projected ellipticity of e ~ 0.4 to 0.5 in the 
inner regions of the galaxy. However, on scales r ~ r v ; r Kazantzidis et al. 1108 ] find that the projected ellipticity is 
£ ~ 0.3. Since both the weak lensing shear and satellite dynamics are determined primarily by the large-scale mass 
distribution of the halos, the roundness of the mass distribution on small scales due to gas cooling should not have 
a dramatic effect. From a theoretical standpoint, therefore, it is not at all unreasonable to expect that galaxy-galaxy 



lensing and satellite dynamics should reflect a significant flattening of the halos. 



Evidence for Flattened Halos from Galaxy-Galaxy Lensing 

Unlike a spherically-symmetric lens for which the gravitational lensing shear is isotropic about the lens center, the 
shear due to an elliptical lens is anisotropic about the lens center. Specifically, at a given angular distance, 9, from 
an elliptical lens, source galaxies which are located closer to the major axis of the mass distribution of the lens will 
experience g reate r shear than sources whic h are located closer to the minor axis (e.g., HI). Noting this, Natarajan 
& Refregier and Brainerd & Wright 1 110] modeled the dark matter halos of field galaxies as infinite singular 
isothermal ellipsoids and made rough estimates of the sizes of observational data sets which would be required to 
detect "anisotropic" galaxy-galaxy lensing and, hence, to constrain the net flattening of the halo population. Both 
studies concluded that, if the mean flattening of the halos is of order 0.3, then only a relatively modest amount of 
imaging data would be necessary to observe the effects of halo flattening on the weak lensing signal. 

In estimating the amount of data that would be required to detect anisotropic galaxy-galaxy lensing, both Natarajan 
& Refregier 1 109] and Brainerd & Wright (H3l 

made the simplifying assumption that each distant source galaxy is 
lensed by only one foreground galaxy. However, for a somewhat deep imaging survey (Iu m ~ 23), the simulations 
of galaxy-galaxy lensing performed by Brainerd, Blandford & Smail 1 13] indicated that most of the galaxies with 
magnitudes in the range 22 < / < 23 would have been lensed at a comparable level by two or more foreground 
galaxies. In a realistic data set, these multiple weak deflections might significantly affect the signal-to-noise that 
could be achieved when attempting to detect anisotropic galaxy-galaxy lensing. This motivated Wright & Brainerd 
ill llll to carry out detailed Monte Carlo simulations of galaxy-galaxy lensing by flattened halos, including the effects 
of multiple weak deflections on the final images of distant galaxies. 

Wright & Brainerd 1 1 1 1 ] showed that multiple weak deflections create systematic effects which could hinder obser- 
vational efforts to use weak lensing to constrain the projected shapes of the dark matter halos of field galaxies. They 
modeled the dark matter halos of lens galaxies as truncated singular isothermal ellipsoids, and for an observational 
data set in which the galaxies had magnitudes in the range 19 < / < 23, they found that multiple deflections resulted 
in strong correlations between the post-lensing image shapes of most foreground-background pairs of galaxies. Im- 
posing a simple redshift cut during the analysis of the data set, zi < 0.5 and z s > 0.5, was sufficient to reduce the 
correlation between the final images of lenses and sources to the point that the expected anisotropy in the weak lens- 
ing signal was detectable via a straightforward average. Wright & Brainerd 1 1 1 1] concluded that previous theoretical 
calculations of weak lensing due to flattened halos had considerably underestimated the sizes of the observational 
data sets which would be required to detect this effect. In particular, for a multi-color survey in which the galaxies 
had apparent magnitudes of 19 < / < 23 and the imaging quality was modest, Wright & Brainerd ll 1 ill found that a 
4<7 detection could be obtained with a survey area of order 22 sq. deg., provided photometric redshift estimates were 
made for the galaxies, the typical error in z p hot was < 0.1, and only source galaxies with azimuthal coordinates that 
were within ±20° of the lens symmetry axes were used in the data analysis. 

To date, only one intrepid team of investigators has claimed a detection of flattened halos from observations of 
galaxy-galaxy lensing. In their analysis of the RCS galaxy-galaxy lensing signal Hoekstra, Yee & Gladders 1 4-8] took 
the approach of modeling the lens galaxies as having halos with ellipticities that scaled linearly with the ellipticity of 
the image of the lens: £h a i = /eiight- Further, they assumed that the major axis of the lens image was aligned with the 
major axis of the halo in projection on the sky. This is a sensible assumption provided the majority of the lenses are 
relaxed systems, and it is justified at least partially by the observations of Kochanek 1 1 12] who found that the major 
axes of the mass and light of strong lens galaxies were aligned to within ~ 10° in projection on the sky. 

Hoekstra, Yee & Gladders |48] performed a maximum likelihood analysis and concluded that spherical halos (i.e., 
/ = 0) could be ruled out at the 99.5% confidence level on the basis of their weak lensing signal (see Figure lot . 
Formally, Hoekstra, Yee & Gladders H found / = 0.77+°;^. Since the mean ellipticity of the lens images in their 
study was (eiight) = 0.414, this implies a mean halo ellipticity of (fihaio) = 0.331qo9 anc ^ a projected axis ratio of 
c/a = 0.67^jj This is in excellent agreement with the expectations for CDM halos, as well as previous observational 
constraints on halo flattening obtained on large physical scales (see, e.g., 1 100]). While it may yet be a bit premature 
to call this result a "definitive" measurement of the flattening of field galaxy halos, it is certainly impressive and the 
statistics will only improve as weak lensing surveys become larger. 




/ 

FIGURE 13. Confidence bounds with which spherical halos can be rejected on the basis of galaxy-galaxy lensing in the RCS 
1 48]. Halos of lens galaxies were modeled as having ellipticities of £h a i = /flight an d the principle axes of the halo mass were 
assumed to be aligned with the symmetry axes of the lens images in projection on the sky. Round halos, / = 0, are excluded at the 
99.5% confidence level. Figure kindly provided by Henk Hoekstra. 

Evidence for Flattened Halos from Satellite Galaxies 

In the case of substantially flattened halos of host galaxies, one would naively expect that satellite galaxies would 
show a somewhat anisotropic distribution about the host. That is, barring possible effects due to infall rates and orbital 
decay, one would expect the satellites to have some preference for being located near to the major axis of the host's 
halo. Until very recently, however, such an observation had not been confidently made and, moreover, a preference for 
clustering of satellite galaxies alon g the m inor axes of host galaxies has been reported at a statistically significant level 
by a handful of authors (| 20], |68], ll 13l0 . The apparent alignment of satellite galaxies with the minor axes of the host 
galaxies is often referred to as the Holmberg effect and in the naive picture of satellite orbits in flattened potentials, 
observations of the Holmberg effect lead to the uncomfortable conclusion that not only is the halo mass flattened, but 
it is also anti-aligned with the luminous regions of the galaxy. 

While o ne is tempted to dismiss the minor axis clustering of satellites observed by Zaritsky et al. l2(ill and 
Holmberg 111 311 as being due to some combination of selection biases and very small sample sizes, it is not easy 
to use this argument for the results of Sales & Lambas |68]. In their study, Sales & Lambas |68] selected hosts and 
satellites from the final data release of the 2dFGRS, with a resulting sample size of 1498 hosts and 3079 satellites. 
The satellites were constrained to be within projected radii r p < 500 kpc of their host and to be within a velocity 
difference \dv\ < 500 km sec -1 . Further, host images were required to have eccentricities of at least 0.1 in order 
that the orientation of their major axes be well-determined. When Sales & Lambas |68] searched their entire sample 
for anisotropies in the distribution of satellites about 2dFGRS hosts, their results were consistent with an isotropic 
distribution. However, when they restricted their sample to only hosts and satellites whose radial velocities differed by 
\dv\ < 160 km sec - 1 , an apparently strong detection of the Holmberg effect (i.e., minor axis clustering of the satellites) 
was found. 

More recently, Brainerd 1 64] investigated the distribution of satellites about hosts in the second data release of the 
SDSS. She selected her samples using three different criteria: (1) the criteria used by Sales & Lambas l68il in their 
investigation of the Holmberg effect for 2dFGRS galaxies, (2) the criteria used by McKay et al. l65il and Brainerd 
& Specian |66] in their analyses of satellite dynamics in the SDSS and 2dFGRS, respectively, and (3) the selection 
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FIGURE 14. Normalized probability distribution of the location of satellite galaxies relative to the major axes of host galaxies in 
the second data release of the SDSS |64]. Dashed line shows the expectation for an isotropic distribution. Formal confidence levels 
at which isotropic distributions can be rejected via x 2 tests are shown in each panel. Also shown is (tj>), the mean value of the angle 
between the major axis of the host galaxy and the direction vector that connects the centroids of the host and satellite. 



criteria used by Zartisky et al. 1 20] in their investigation of the Holmberg effect. In addition, Brainerd 1 64] restricted 
the analyses to hosts with ellipticities e > 0.2 and satellites that were found within a projected radius of 500 kpc. The 
three selection criteria lead to samples of: (1) 1351 hosts and 2084 satellites, (2) 948 hosts and 1294 satellites, and (3) 
400 hosts and 658 satellites respectively. 

In all three samples, Brainerd |64| found that the distribution of satellites about their hosts was inconsistent with an 
isotropic distribution. Formally, when a Kolmogorov-Smirnov test was applied to the distribution of satellite locations, 
an isotropic distribution was rejected at a confidence level of > 99.99% for sample 1, > 99.99% for sample 2, and 
99.89% for sample 3. Further, the mean angle between the major axis of the host and the direction vector on the sky that 
connected the centroids of the hosts and satellites was found to be (<j>) = 41. 6° ±0.6° for sample 1, (<j>) =41.6°±0.7° 
for sample 2, and (<p) = 41 .6° ± 1 .0° for sample 3. That is, a clear anisotropy in the distribution of satellites about the 
hosts was seen, and the satellites showed a preference for being aligned with the major axis of the host rather than the 
minor axis (see Figure \14\. In addition, Brainerd 1 64] investigated the dependence of (<j>) with projected radius on the 
sky and found that the majority of the anisotropy arose on small scales (< 200 kpc) in all three samples (see FigurefOJ. 
In other words, the anisotropy was detected on physical scales that are comparable to the expected virial radii of large, 
bright galaxies. On scales much larger than the expected virial radii of galaxy-sized halos (r p ~ 400 kpc to 500 kpc), 
the distribution of satellites about the SDSS hosts was consistent with an isotropic distribution at the la level. 

Aside from the Brainerd 16411 claim of "planar" (rather than "polar") alignment of satellites with the symmetry axes 
of their hosts, there has been only one other similar claim. Valtonen, Teerikorpi & Argue ll 14ll found a tendency for 
compact satellites to be aligned with the major axes of highly-inclined disk galaxies; however, their sample consisted 
of only 7 host galaxies. Although it is extremely tempting to accept its veracity based upon an intuitive sense that 
planar alignment of satellites is more dynamically sensible than polar alignment, it is clear that the Brainerd |64| 
result is badly in need of independent confirmation. 

Sales & Lambas |68] used a data set of very similar size (and in one case identical selection criteria) to that of 
Brainerd |64| yet did not detect any anisotropy in the satellite distribution when they analysed their entire sample. 
Why this is the case remains a mystery at the moment, but it may be attributable to a combination of two things. First, 
the velocity errors in the 2dFGRS are typically larger than those in the SDSS (~ 85 km sec _1 versus ~ 20 km sec~'to 
~ 30 km sec ). At some level, this would lead to a higher fraction of interlopers (i.e., false satellites) in the Sales 
& Lambas |68] sample than in the Brainerd 1 64] samples. Second, van den Bosch et al. ll 1 3fl found that when they 
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FIGURE 15. Mean orientation of satellite galaxies with respect to the major axes of the hosts as a function of the projected 
radius for galaxies in the second data release of the SDSS 1 64 ] . Dashed line shows the expectation for an isotropic distribution. 
Here h = 0.7 has been adopted. 



combined mock redshift surveys with the 2dFGRS, there was a clear absence of satellites at small projected radii in 
the 2dFGRS. Since the majority of the anisotropy seen by Brainerd 1 64] appears to come primarily from small scales, 
it could be that Sales & Lambas l68ll simply had too few pairs of hosts and satellites at small separations to detect the 
anisotropy. Any lack of host-satellite pairs in the 2dFGRS data, however, does not explain why a Holmberg effect was 
detected by Sales & Lambas |68] when they restricted their analysis to host-satellite pairs with \dv\ < 160 km sec . 
When Brainerd 1 64] imposed the same restriction on her sample 1 (i.e., the sample selected using the Sales & Lambas 
1 68 ] selection criteria), she found that the satellites with \dv\ < 160 km sec -1 displayed an anisotropy that was identical 
to that of the full sample: a clear alignment of the satellites with the host major axes. The cause of this discrepancy 
is not at all obvious. It may in part be attributable to the fact that a value of \dv\ = 160 km sec~'is comparable to the 
error in a typical measurement of \dv\ for hosts and satellites in the 2dFGRS. Also, work by van den Bosch et al. 1 69] 
suggests that the interloper fraction is substantially higher for host-satellite pairs with low values of \dv\ than it is for 
host-satellite pairs with high values of \dv\. It could, therefore, be possible that the Sales & Lambas 1 68 ] sample with 
\dv\ < 160 km sec~'is heavily contaminated with interlopers and some strange, unknown selection bias is giving rise 
to their signal. 

Finally, it is worth notating that not only are the observational conclusions about the distribution of satellite galaxies 
particularly muddy at the moment, so too are the theoretical conclusions. Zaritsky et al. 1 20] compared their observed 
Holmberg effect with high-resolution CDM simulations and were unable to recover their observations. Penarrubia et 
al. II 1611 investigated both polar and planar orbits of satellites inside a massive, flattened dark matter halo and found 
that the planar orbits decayed more quickly that the polar orbits. They therefore suggest that such differences in orbital 
decay rates could be the origin of the Holmberg effect. Abadi et al. II 1711 suggest that the Holmberg effect could be 
caused by the cumulative effects of accretion of satellites by the primary. However, Knebe et al. 1 1 18] found that the 
orbits of satellites of primary galaxies in cluster environments were located preferentially within a cone of opening 
angle 40° (i.e., planar alignment, n ot po l ar). S ince the structure of cold dark matter halos is essentially independent 
of the mass scale of the halo (e.g., ill 1 9fl . Il20l0 . the implication of this result would be a preference for the satellites 
of isolated galaxies to be aligned with the major axis of the host. All of this in mind, perhaps the only answer to 
the question 'Are either the Sales & Lambas 1 68 ] or Brainerd 1 64 ] observations of anisotropic satellite distributions 
consistent with galaxy halos in a CDM universe?" is, for now, "Maybe". 



SUMMARY 



There has been a long period of time over which it has been perfectly acceptable to write papers on investigations into 
the nature of the dark matter halos of field galaxies that begin with a statement along the lines of "Although modern 
theories of galaxy formation posit that all large galaxies reside within massive halos of dark matter, the characteristic 
properties of those halos (e.g., mass, radial extent, and shape) are not well-constrained by the current observations". 
That time is now coming to an end. The wealth of data that has been acquired in recent years is truly beginning to 
place strong, direct constraints on the dark matter halos of field galaxies. 

Weak lensing and satellite dynamics have proven themselves to be excellent probes of the gravitational potentials 
of large, bright galaxies on physical scales r > 100 h~ kpc. While one might be skeptical and discount the results 
that come from one technique or the other, the fact that both are yielding consistent constraints cannot be ignored. 
Both weak lensing and satellite dynamics lead to statistical constraints on the halo population as a whole, rather than 
constraints on any one particular galaxy halo, and it is especially the acquisition of extremely large data sets that 
has allowed these techniques to begin to fulfill their promise of mapping out the gravitational potentials associated 
with large, massive halos. Weak lensing and satellite dynamics have inherent advantages and disadvantages, but since 
their systematic errors and selection biases are completely uncorrected, they are extremely complementary to each 
other. At least at the moment, when strong constraints are only just beginning to emerge from each technique, this 
complementarity is very reassuring. 

Based upon my own critical, and hopefully unbiased, reading of the recent literature, I think it is fair to say that, 
both individually and in combination, weak lensing and satellite dynamics are pointing toward the following scenario 
for the nature of large, bright field galaxies and their halos: 

• The dark matter halos are well-characterized by NFW-type objects in terms of their gravitational properties. The 
dynamics of satellite galaxies strongly prefer NFW halos to isothermal halos. 

• The virial masses that are inferred for large field galaxies are in good agreement with the predictions for galaxy- 
mass halos in the context of cold dark matter. Specifically, the virial mass of the halo of an "average" L* galaxy 
is in the range (8 — 10) x lO n h~ 1 M <3 when NFW profiles are fit to the data. 

• There are clear differences in the depths of the potential wells of the halos that surround galaxies of differing 
morphology and differing intrinsic luminosity. Specifically, the virial masses of the halos of L* ellipticals exceed 
those of L* spirals by a factor of at least 2. The actual value of the mass excess depends upon details of the 
data and its analysis. In addition, the virial masses of the halos of high luminosity galaxies exceed those of low 
luminosity galaxies. Again, however, the amount by which they differ depends upon details of the data and its 
analysis. 

• Averaged over all galaxies with L>L* , the mass-to-light ratio computed on scales larger than the optical radii of 
the galaxies is, at most, weakly-dependent upon the luminosity of the galaxy. At the 2a level, the mass-to-light 
ratio of the average galaxy with L>L* is consistent with a constant value. 

• The dark matter halos are flattened, rather than spherical, and the degree of flattening on large scales (<~ 100 kpc 
to <~ 200 kpc) is consistent with the predictions of cold dark matter. 

It is worth noting that the above list comes from quite a diverse set of data. In particular, the data are spread over a 
wide range in redshift. With the exception of preliminary data from DEEP2, the satellite dynamics studies have median 
redshifts of z mec i <~ 0.07. The weak lenses in the SDSS data have a median redshift of z me d ~ 0.16 and the weak lenses 
in the RCS and COMBO-17 data have considerably higher redshifts, z me( j ~ 0.4. Since it is clear that the field galaxy 
population has evolved since z ~ 0.5, it is not entirely fair to lump the results from all of these studies together, and 
I think the big challenge to the weak lensing community in particular will be to eventually place constraints on the 
evolution of field galaxies and their halos from, say, z ~ 1 to the present. 

Nevertheless, I think we have reached a particularly gratifying time in which we are really being able to measure 
some of the fundamental properties of dark mater halos on physical scales that extend well beyond the visible images 
of the galaxies at their centers. A remarkably consistent picture of the large-scale gravitational properties of the halos 
is emerging from the observations and, at least for now, that picture seems entirely in accord with a cold dark matter 
universe. 
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